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Using 230.2 fb _1 of e + e~ annihilation data collected with the BABAR detector at and near the peak 
of the T(45) resonance, 489 ± 55 events containing the pure leptonic decay Df — » have been 

isolated in charm-tagged events. The ratio of partial widths T(Df —> ^ i/ M ) /T(Df — > 4>tt + ) is 
measured to be 0.143 ± 0.018 ± 0.006 allowing a determination of the pseudoscalar decay constant 
fo s = (283 ± 17 ± 7 ± 14) MeV. The errors are statistical, systematic, and from the Df — » 4>tt + 
branching ratio, respectively. 

PACS numbers: 13.20.He, 14.40.Nd, 14.60.Fg 



Measurements of pure leptonic decays of charmed pseu- 
doscalar mesons are of particular theoretical importance. 
They provide an unambiguous determination of the over- 
lap of the wavefunctions of the heavy and light quarks 
within the meson, represented by a single decay constant 
(/m) for each meson species (M). The partial width for 
a £)+ meson to decay to a single lepton flavor (I) and its 
accompanying neutrino (yi), is given by 

where m D and m ; are the Df and lepton masses, re- 
spectively, Gf is the Fermi constant, and V cs is the 
CKM matrix element giving the coupling of the weak 
charged current to the c and s quarks The partial 
width is governed by two opposing terms in mf. The 
first term reflects helicity suppression in the decay of the 
spin-0 meson, requiring the charged lepton to be in its 
unfavored helicity state. The second term is a phase- 
space factor. As a result, the ratio of r : fi : e de- 
cays is approximately 10 : 1 : 0.00002. Lattice calcula- 
tions have resulted in fp s = (249 ± 17) MeV and a ratio 
f Ds /f D = 1.24±0.07 2]. CLEO-c has recently measured 
a value for f D = (223 ± 17) MeV 0]. 

We present herein the most precise measurement to 
date of the ratio r(L>+ /j,+u^)/T(D+ -> (jm + ) and 
the decay constant fjj 3 . The data (230.2 fb -1 ) were col- 
lected with the BABAR detector at the asymmetric-energy 
e + e~ storage ring PEP-II at and below the T(45) reso- 
nance. The BABAR detector is described in detail else- 
where 0. Briefly, the components used in this analysis 
are the tracking system composed of a five-layer silicon 
vertex detector and a 40-layer drift chamber (DCH), the 
Cherenkov detector (DIRC) for charged tt— K discrim- 
ination, the CsI(Tl) calorimeter (EMC) for photon and 
electron identification, and the 18-layer flux return (IFR) 
located outside the 1.5 T solenoid coil and instrumented 
with resistive plate chambers for muon identification and 
hadron rejection. 

The analysis proceeds as follows. In order to measure 
D+ — > /x+i^, the decay chain D* s + — * 7_D+,D+ — * 
is reconstructed from D* + mesons produced in the hard 



fragmentation of continuum cc events. The subsequent 
decay results in a photon, a high-momentum and 
daughter muon and neutrino, lying mostly in the same 
hemisphere of the event. Signal candidates are required 
to lie in the recoil of a fully reconstructed D°, D + , Df , or 
D* + meson (the "tag"), wherein the tag flavor, and hence 
the expected charge of the signal muon, is uniquely deter- 
mined. To eliminate signal from B decays, the minimum 
tag momentum is chosen to be close to the kinematic 
limit for charm mesons arising from B decays. Tagging 
in this manner significantly reduces backgrounds, while 
improving the missing mass resolution of the signal. 

Tag candidates are reconstructed in the follow- 
ing modes: D° — » K ~it + (it ), K~ir + Tr + ir~ , D + — > 
A'-7r+7r+(7r°), K°n + (ir ), K°-k+ tt+ tt~ , K+K~tv + , K° S K+ , 
D+ -> K° s K+. / cj)p + , and D*+ -> n+D°, with D 6 -> 
A!?7r+7r-(7r ), K°K+K~, A°7r°. Kaons are identified 
using information from the DCH and the DIRC. Re- 
quirements on the vertex probability of the tag decay 
products are imposed. For each tag mode a signal region 
and sideband regions in the tag mass distribution are 
defined. The signal region spans ±2 standard deviations 
(ctag) around the mean (/itag), determined from fits to 
the tag mass distribution in data events. The sidebands 
extend from 3 to 6 er tag on either side of /x t ag (Fig-D- 

For each event a single tag candidate is chosen and then 
used in the subsequent analysis. To pick this tag among 
multiple candidates within an event (there are 1.2 candi- 
dates on average in events with at least one candidate) 
modes of higher purity are preferred. In events where 
two tag candidates are reconstructed in the same mode, 
the quality of the vertex fit of the D meson is used as 
a secondary criterion. After subtracting combinatorial 
background there are 5 * 10 5 charm tagged events with a 
muon amongst the recoiling particles. 

The signature of the decay D* + — > jD+ is a narrow 
peak in the distribution of the mass difference AM = 
M(iivi)-M(p,v) at 143.5 MeV/c 2 . The D*+ signal is re- 
constructed from a muon and a photon candidate in the 
recoil of the tag. Muons are identified as non-showering 
tracks penetrating the IFR. The muon must have a mo- 
mentum of at least 1.2 GeV/c in the center-of-mass (CM) 
frame and have a charge consistent with the tag flavor. 
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FIG. 1: Tag mass distribution, showing the signal and side- 
band regions, in events with a recoil muon. All tag modes 
are combined, scaling their mass and width to that of the 
D° -> R-TT+ mode. 



Muons used in this analysis are identified with an average 
efficiency of sa 70 %, while the pion misidentification rate 
is « 2.5 %. Clusters of energy in the EMC not associated 
with charged tracks are identified as photon candidates. 
The photon CM energy must exceed 0.115 GeV. 

The CM missing energy (-B^lss) anc ^ momentum (p^ liss ) 
are calculated from the four-momenta of the incoming 
e + e~, the tag four-momentum, and the four-momenta 
of all remaining tracks and photons in the event. The 
energy of the charged particles that do not belong to 
the tag is calculated from the track momentum under 
a pion mass hypothesis. Assigning a mass according to 
the most likely particle hypothesis has negligible effect 
on the missing energy resolution. Since the neutrino in 
the signal decay leads to a large missing energy in the 
event, the requirement E^ iss > 0.38 GeV is made. 

The neutrino CM four- momentum (p* = (\Pt\,pt)) 
is estimated from the muon CM four-momentum (p*) 
and p^ iss , using a technique adopted from Ref. |j§. The 
difference |pJnj ss — p^. | is minimized, while the invariant 
mass of the neutrino-muon pair is required to be the 
known mass of the Df |(|. Studies of simulated decays 
of signal and background cc events show that the quan- 
tity p C orr = |p?ni S sl — \pt\ 1S centered at for signal de- 
cays, while for the cc events it peaks at a negative value 
significantly separated from the signal. A requirement 
Peon- > —0.06 GeV/c is imposed. To reduce contributions 
from background events where particles are lost along the 
beam pipe in the forward direction, a requirement on the 
neutrino CM polar angle 9* > 38 is made. The muon 
CM four-momentum (p*) is combined with p* to form 
the Df candidate. Unlike the signal D+, a large num- 
ber of random Df combinations have the muon candi- 
date aligned with the Df flight direction. A requirement 
cos(a fJi ,D, ) < 0.90 is made on the angle between the muon 
direction in the Df frame and the flight direction in 
the CM frame. The Df candidate is then combined with 
a photon candidate to form the -D* + . The CM momen- 
tum of correctly reconstructed D* + is typically higher 



than that of random combinations; signal candidates are 
required to have > 3.55 GeV/c. The resulting sig- 

nal detection efficiency in tagged events is esi g = 8.13%. 

The selection requirements on E^ iss , Oin,D„, Pcorr, 0*, 
and |j9^„+ 1 are optimized using simulation to maximize 

the significance s/y s + b, where s and b are the sig- 
nal and background yields expected in the data set. 
Backgrounds arise from several distinct sources. The 
first class of background are events e + e~ — > //, where 
/ = u,d, s, 6, or t, which do not contain a real charm tag. 
The contribution of these events is estimated from data 
using the tag sidebands. In addition there are events 
e + e~ — > cc where the tag is incorrectly reconstructed. 
Although these events potentially contain the signal de- 
cay, they are also subtracted using the tag sidebands. 
These two sources amount to w 42 % of the background. 

The second class of background events (« 26 %) are 
correctly tagged cc events with the recoil muon coming 
from a semileptonic charm decay or from t + — ► [i^v^Vr. 
This includes events D* + — > jDf — ► 7T + ^ T , t + — > 
H+v^Dt. To estimate the size and shape of this back- 
ground contribution, the analysis is repeated, substitut- 
ing a well-identified electron for the muon. Except for a 
small phase-space correction, the widths of weak charm 
decays into muons and electrons are assumed to be equal. 
QED effects such as bremsstrahlung (e + — > 7e + ) energy 
losses and photon conversion (7 — > e + e~), where the 
muon equivalents have a much lower rate, are explicitly 
removed. In particular, bremsstrahlung photons found 
in the vicinity of an electron track are combined with 
the track. The small number of events with an elec- 
tron from a converted photon that survive the selection 
are suppressed by a photon conversion veto, using the 
vertex and the known radial distribution of the mate- 
rial in the detector. The muon selection efficiency as a 
function of momentum and direction is measured using 
e + e~ — * /i + /i _ 7 events, while radiative Bhabha events 
are used to quantify the electron efficiency. The ratio 
of muon to electron efficiencies is applied as a weight to 
each electron event. 

The remaining backgrounds are estimated from simu- 
lation. These include events (« 20 %) with pure leptonic 



where the 



decays of a Df or D + meson, D 

Df g s is produced either directly in cc fragmentation or in 
decays of D1\, excluding the signal decay chain. If the 
photon used in the reconstruction originates from a 7r° of 
a D*, \ decay, the AM distribution peaks sharply around 



W 

70 MeV/c 2 ; otherwise it is flat, 
(w 1 %) arises from decays D* + 



A small background 
> "/Df — > r yT + v T with 
(tt°)v t and the charged pion being misidcntificd 
as a muon. Its AM distribution peaks close to that of 
the signal. Other backgrounds (f» 10 %) include signal 
events with an incorrectly chosen photon candidate, and 
hadronic cc events with one of the final state hadrons, 
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FIG. 2: AM distribution of charm-tagged events passing the 
signal selection. The tag can be from the tag signal region 
(solid lines) or the sidebands (dashed lines). In the bottom 
plot the signal muon is replaced with an electron to estimate 
the semileptonic charm and r decay background. 



FIG. 3: AM distribution after the tag sidebands and the 
electron sample are subtracted. The solid line is the fitted 
signal and background distribution (iVsig/sig + iVskgd /Bkgd ) , 
the dashed line is the background distribution (A^Bkgd/Bkgd) 
alone. 



usually a tt + or a K + , being misidentified as a muon. 
These backgrounds have a flat AM distribution. 

Events that pass the signal selection are grouped into 
four sets, depending on whether the tag lies in the sig- 
nal region or the sideband regions, and on whether the 
lepton is a muon or an electron (Fig. [5J . For each lepton 
type the sideband AM distribution is subtracted. The 
electron distribution, scaled by the relative phase-space 
factor (0.97) appropriate to semileptonic charm meson 
decays and leptonic r decays is then subtracted from 
the muon distribution. The resulting AM distribution 
is fitted with a function (iV S ig/sig + A r Bk g d/Bk g d)(AM), 
where /si g and /Bkgd describe the simulated signal and 
background AM distributions. The function /gig is a 
double Gaussian distribution. The function /Bkgd con- 
sists of a double and a single Gaussian distribution de- 
scribing the two peaking background components, and 
a function describing the flat background component. 
The relative sizes of the background components, along 
with all parameters except iVsig and A^kgd are fixed to 
the values estimated from simulation. The x 2 fit yields 
A'sig = 489 ± 55 (stat) signal events and has a fit proba- 
bility of 8.9% (Fig. El . 

The branching fraction of Df — > ^ + v^ cannot be de- 

(*)+ 

tcrmincd directly, since the production rate of D s 
mesons in cc fragmentation is unknown. Instead the par- 
tial width ratio r(Z3+ — > n + v^) /T(D+ — ► <fiir + ) is mea- 
sured by reconstructing D* + — > 7-D+ — > 707r + decays. 
The Df — > [i^v^ branching fraction is evaluated using 
the measured branching fraction for D+ — > 4>ir + . 

Candidate <j) mesons are reconstructed from two kaons 
of opposite charge. The <j> candidates are combined with 
charged pions to form Z)+ meson candidates. Both times 
a geometrically constrained fit is employed, and a mini- 
mum requirement on the fit quality is made. The 4> and 
the -D+ candidate masses must lie within 2 a of their 
nominal values, obtained from fits to simulated events 
and data. Photon candidates are then combined with 



the Df to form D* + candidates. The same require- 
ments on the CM photon energy and D* + momentum 
as in the — > signal selection are made. The 

D* + — > -fDf — > selection efficiency in tagged 

events is = 9.90%. Data events that pass the selec- 
tion are grouped into two sets: the tag signal and side- 
band regions. After the tag sideband has been subtracted 
from the tag signal AM distribution, the remaining dis- 
tribution is fitted with (iV^/^+iY^Bkgd/^BkgdXAM), 
where is a triple Gaussian, describing the simu- 
lated D* + — ► 7D+ — > j4>tt + signal, and /^Bkgd con- 
sists of a broad Gaussian centered at 70 MeV/c 2 and 
a function Q describing the simulated background AM 
distributions. The Gaussian describes the background 
D* + — > TT°Df — > 7r°07r + where the photon candidate 
originates from the ir°. The relative sizes of the back- 
ground components, along with all parameters except 
Nfa, -/V^Bkgdj and the mean of the peak are fixed to 
the values estimated from simulation. The x 2 fit yields 
Ntfrn — 2093 ± 99 events and has a probability of 25.0% 
(Fig. El . From simulation 48 ±23 events D*+ — > jD+ — ► 
r yfo(980)(K + K^)n + are expected to contribute to the 
signal, where the error is mostly from the uncertainty in 
the D+ -> fo(980)(K+K-)Tr+ braching ratio. 

Precise knowledge of the efficiency of reconstructing 
the tag is not important, since it mostly cancels in the 
calculation of the partial width ratio. However, the pres- 
ence of two charged kaons in Df — ► <fiir + events leads 
to an increased number of random tag candidates, com- 
pared to D+ — ► events, which decreases the chances 
that the correct tag is picked. The size of the correction 
for this effect to the efficiency ratio (e^/esig) is deter- 
mined to be —1.4% in simulated events. 

To measure the effect of a difference between the 
D* + momentum spectrum in simulated and data events, 
D* + — ► "iD+ — > 707r + events are selected in data with 
the D* + momentum requirement removed. The sample is 
purified by requiring the CM momentum of the charged 
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FIG. 4: AM distribution of selected D* + -> 7D+ -> 7^71-+ 
events after the tag sideband is subtracted. The solid line 
is the fitted signal and background distribution (Njm + 
A^TrBkgd/^vrBkgd), the dashed line is the background distribu- 
tion (JV^Bkgd/^TrBkgd) alone. 



pion to be at least 0.8 GeV/c. The efficiency-corrected 
D* + momentum distribution in data is compared to that 
of D* + in simulated D* + — ► jDf — > 707T+ events. A 
harder momentum spectrum is observed in data. The de- 
tection efficiencies for signal and £>* + — > jDf — > j<pn + 
events are re-evaluated after weighting simulated events 
to match the D* + momentum distribution measured in 
data. The correction to the efficiency ratio is +1.5%. 

With both corrections applied, the partial width ra- 
tio is determined to be T^/T^ — (N/e)si g /(N/e),j )7T x 
B(4> -> K+K-) = 0.143 ± 0M8(stat), with B(4> -> 
K+K-) = 49.1% 0. 

The combined systematic uncertainty due to the cor- 
rections applied, taken as half the size of each correction, 
is 1.0%. The systematic error in the signal efficiency due 
to selection criteria insensitive to the D* + momentum 
is evaluated using reconstructed D*° — > -fD° — > jK~n + 
events. The conditions present in the signal are emulated 
by removing the charged pion, taken to represent the neu- 
trino, from these events. The signal reconstruction and 
selection steps are repeated, and the selection efficiencies 
compared between simulated and data events. The as- 
signed systematic uncertainty is 1.4%. For the D+ — > 
4>tt + selection, requirements on the Dj and <fi vertex fit 
probability contribute a systematic uncertainty of 0.7%, 
estimated from comparisons of Df — > </>7r + events in sim- 
ulation and data. Control samples of e + e~ — > /j, + fi~j 
and D* + — ► ir + D° — ► n + K~ir + events are used to mea- 
sure the particle identification efficiencies of muons and 
charged kaons and pions in data, and to correct the sim- 
ulated signal and D* + — + 'yDf — + ^c/)tt + efficiencies. An 
uncertainty of 0.7% is associated with these corrections, 
mainly due to the limited statistics of the control sam- 
ples. The systematic uncertainties in the track recon- 
struction efficiency cancel partially in the Df — > fi Vp, 
to — > <f)ir + ratio and contribute 1.2%. An additional 
uncertainty of 1.1% is due to the statistical limitations 
of the simulated signal and Df — > (fm + event samples. 



Simulation studies are used to evaluate the systematic 
uncertainties arising from a possible inadequate param- 
eterization of the signal (0.9 %) and background (2.3 %) 
shapes. Simulations are also used to determine the sys- 
tematic uncertainty associated with the subtraction of 
the electron sample (0.4%). The error on the branching 
ratio B(4> ~ * K + K~) is 1.2%, the uncertainty on the 
D+ -> /o(980)7r+ background is 1.1%. The total sys- 
tematic uncertainty on T(Df — > ^, + u fl )/T(Df — > (f)ir + ) 
is 3.9%. 

Using the BABAR average for the branching ratio 
B(D+ -> 4>tt+) = (4.71 ± 0.46)% we obtain the 

branching fraction B(D+ -> fi+v^) = (6.74±0.83±0.26± 
0.66) x 10~ 3 and the decay constant f D , = (283 ± 17 ± 
7 ± 14) MeV. The first and second errors are statistical 
and systematic, respectively; the third is the uncertainty 
from B(D+ — > 4>ir + ). The ratio of our value for fo s to fo 
from the CLEO-c measurement, JdJ/d = 1.27 ± 0.14, 
is consistent with lattice QCD. 

Using B{D+ -► 07r+) PDG = (3.6 ± 0.9) % 0, the 
branching fraction is B(D+ — ► /i + ^) = (5.15 ± 0.63 ± 
0.20 ± 1.29) x 10~ 3 and the decay constant f Ds = (248 ± 
15 ± 6 ± 31) MeV. 
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